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EARTH-MOON SYSTEM: DYNAMICS AND PARAMETER ESTIMATION 

By 

W. J. Breedlove, Jr. 1 

INTRODUCTION AND SUMMARY 
General 

This report constitutes the final report required under 
grant NSG 1152 issued by the National Aeronautics and Space 
Administration entitled "Earth-Moon System: Dynamics and Param- 

eter Estimation." A preliminary survey was prepared under 
Master Contract Agreement NAS1-11707, Task Authorization 
No. 41. The results of that survey are in reference 1. 

Original Objectives 

Objectives of the original grant were to 

1. Prepare a dynamic model of the coupled translational/ 
rotational motion of the moon for use in a data-reductxon 
process at NASA/Langley Research Center (LaRC) . The_ available 
Lunar Laser Ranging Experiment (LURE) data perhaps in conjunc- 
tion with other data types were to be utilized; 

2. Obtain an accurate numerical integration scheme 
capable of solving the equations of motion to a high degree 
of accuracy; 

3. Develop observational and normal equations for the 
LURE data type; and 

4 . Participate in the data analysis and interface 
grant-developed software with existing LaRC software for 
use in a parameter-estimation process. 

Associate Professor, Department of Mechanical Engineering 
and Mechanics, Old Dominion University, Norfolk, Virginia, 
23508. 
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Objectives of Continuation Grant 

Objectives in an extension to the original grant included 

5. Conversion of grant-developed software from the 

CDC 6600/6700 computers to the Old Dominion University DEC 10 
coiriputer ; 

6 . Incorporation of BIH data into the software for Earth 
rotation; 

7 . Analysis of LURE data to determine or study the 
possibility of determining various parameters; and 

8. Study of potentialities of using LURE data with other 
data types for parameter estimation. 

Accomplishments 

The following provides an assessment of the status of 
work towards the above objectives. 

1. The Unified Model of Lunar Translation/Rotation 
(UMLTR) has been coded and verified. The model is capable 
of integrating all planets in a coupled sense with the 
rotational motion of the Moon. Comparisons of the UMLTR 
with DE-96 show accuracies of « 500 m for the Earth-Moon 
barycenter after 3 years and w 100 m for the heliocentric 
position of the Moon for the same time span. Comparisons 

of the UMLTR with Eckhardt's 300 series for lunar librations 
show standard deviations of less than 0.15 seconds m the 
physical librations p, la, and t. 

2. A paper "A Unified Model of Lunar Translation and 
Rotation for Use in the Reduction of Lure Data" was presented 
to the Division of Dynamical Astronomy of the AAS . 

3. An accurate integrator was obtained, verified, and 
is currently in use in the UMLTR. The method is "Everhart’s 
Single Sequence Method" (ref. 2) . 

4. A paper, "A Numerical Study of the Effects of Fourth 
Degree Terms in the Earth-Moon Potential on Lunar Physical 
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Librations , " was presented at the Scientific Applications of 
Lunar Laser Ranging Symposium and published in the proceedings 
(Appendix A) . 

5. The BIH data has been obtained from the Jet Propul- 
sion Lab (JPL) ; and has been incorporated in the software. 

6. A program, N0RM, has been obtained from NASA/LaRC, 

This program forms normal equations to be used in the data- 
reduction process. It essentially solves the light-time 
problem, calculates the observable quantity, and forms observa- 
tional related partial derivatives; subroutine T0P0S and 
certain time conversion routines from N0RM have been verified. 

7 . Programming of the variational equations is about 
75 percent complete. Numerical partial derivatives were 
utilized for the entire study. 

8. LURE normal points (and filtered data) for a six- 
year period are available (ref. 3) . 

9. A paper "A Unified Model of Lunar Translation/Ro- 
tation" has been submitted to the Celestial Mechanics Jour- 
nal (Appendix B) . 

10. A study of the coupling between lunar translational 
and rotational motions was performed using the UMLTR. It was 
found that fourth degree lunar harmonics produced orbital 
perturbations of 5 cm, 0.0018 cm, and 0.0025 cm, respectively 
in semimajor axis, perigee, and longitude after three years 
due to the coupling mechanism. 

11. A paper "On Lunar Orbital-Rotational Coupling and 
Figure-Figure Interactions in the Earth-Moon System" has 
been accepted for the 17th Aerospace Sciences Meeting of 
the AIAA (Appendix C) . 

12. A study of the feasibility of estimating simulta- 
neously the lunar inertias A/MR 2 , B/MR 2 , and C/MR 2 using 
actual LURE observation times (ref. 3,) , and numerical partials 
to form a normal matrix has been accomplished. The results of 
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this study show no major problem in estimating these quantities 
when a rather limited set of parameters is considered. The 
results of this study are presented in Appendix C. 

THE UNIFIED MODEL OF LUNAR TRANSLATION/ROTATION (UMLTR) 

A numerical model which integrates the heliocentric 
equations of motion for all planets, the Earth-Moon bary- 
center, the geocentric motion of the Moon, and the Moon's 
rotational motion in a coupled sense has been derived, coded, 
and verified. A description of the UMLTR with details on 
various comparisons with Eckhardt's analytical model for the 
lunar physical librations and the JPL DE-96 ephemeris is 
presented in Appendix B. Reference 4 contains the derivations 
of the basic equations used in the UMLTR. A description of 
the numerical integration package used in the UMLTR as well as 
a description of an auxilliary program ANEAM0 is provided in 
reference 5. Program ANEAM0 evaluates analytic results for 
the rotational and translational motion of the Moon and inter- 
polates DE-96 for the position of the Moon and planets. 
Reference 6 provides a description of some modifications to 
the UMLTR. Also, that reference provides a data analysis 
framework and the complete set of variational equations. A 
description of the various data files, programs, and auxilliary 
programs necessary for use with the UMLTR is presented in 
reference 7. Reference 8 describes the normal equation formu- 
lation to be used with the UMLTR. Finally, the general status 
of the work on NASA Grant NSG-1152 is presented in reference 9. 

THE EFFECT OF FIGURE -FIGURE INTERACTIONS 
ON LUNAR PHYSICAL LIBRATIONS 

One of the main purposes for development of the UMLTR 
was to investigate certain small effects on lunar rotation 
that might be discernible in the LURE data. One such effect, 
the figure-figure interaction (FFI) between the Earth and the 
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Moon, stems from certain terms in the Earth-Moon mutual poten- 
tial which have been neglected in most models to date. 

A disadvantage of numerical methods is that long-term 
behavior cannot be evaluated due to roundoff error growth. 

Thus, the effect of FEI was made only for a three-year period. 
Some preliminary numerical results are presented in Appendix 
A. There it was concluded that the effect of FFI on lunar 
physical librations amounted to at least 0.01" which mapped 
into range errors of 3 to 4 cm. This effect should be modeled, 
therefore, since the LURE data accuracy is expected to approach 
the 2- to 3-cm level. 

Yoder in reference 10 subsequently showed, using an 
approximate theoretical analysis , that the effect of the FFI 
was to produce an 18.6-year libration of magnitude 0.0018" 
in longitude and 0.08" in latitude. This effect is now 
being included in the JPL physical l-ibration model. 

THE EFFECT OF TRANSLATIONAL/ROTATIONAL COUPLING 
ON THE LUNAR ORBIT 

Another small effect heretofore neglected in analytical 
and numerical models is the rotational/translational coupling 
of lunar rotation with the lunar orbit. This effect is 
usually handled in an iterative manner, i.e., an ephemeris is 
produced assuming Cassini’s laws are valid or that a given 
analytical theory represents the lunar rotation. Then, a 
rotational integration is made providing better knowledge for 
the next ephemeris generation. 

The coupling effect is difficult to evaluate numerically. 
To gain some insight into the nature of, and magnitude of, 
the coupling effects, some numerical experiments were made. 
These experiments are described in Appendix C. 

It was found that fourth degree lunar harmonics produce 
perturbations of 5 cm, 0.0018 cm and 0.0025 cm respectively 
in semimajor axis, perigee and longitude, after three years, 
due to the coupling mechanism. 
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AN ERROR ANALYSIS FOR ESTIMATING LUNAR 
INERTIAS FROM LURE DATA 


The feasibility of estimating lunar inertias A/MR 2 , 

B/MR 2 , and C/MR 2 using the UMLTR and a single data set was 
studied. The results of this study are presented in Appendix 
C. This study did not involve a complete parameter set but 
indicates no major correlations or problems in estimating 
the above parameters using a data set such as LURE. 



APPENDIX A 


A NUMERICAL STUDY OF THE EFFECTS OF FOURTH DEGREE 
TERMS IN THE EARTH-MOON MUTUAL POTENTIAL ON 
LUNAR PHYSICAL LIBRATIONS 


(Article from Scientific Applications of Lunar Laser Ranging . ) 


7 



A NUMERICAL STUDY OF THE EFFECTS OF FOURTH DEGREE TERMS IN THE 
EARTH-MOON MUTUAL POTENTIAL ON LUNAR PHYSICAL LIBRATIONS* 


W.J. Breedlove, Jr. 

Mechanical Engineering and Mechanics Department 
Old Dominion University 
Norfolk, Virginia, U.S.A. a 


ABSTRACT. Lunar laser Ranging Data (LURE) is currently available 
with a range accuracy of 10 cm. Accuracies of 2 to 3 cm are 
ultimately expected. Previous analyses have shown that for an 
Earth-Moon range precision of 10 cm, certain fourth degree lunar 
harmonic gravity coefficients must be included in the lunar 
physical libration model. An order-of -magnitude calculation 
shows that some torques arising from previously neglected fourth- 
degree "figure- figure interaction" terms of the Earth-Moon mutual 
potential could produce effects on the physical librations that 
should be modeled for range accuracies of 2 to 3 cm. 

An investigation of the effect of these "interaction" terms has 
been undertaken utilizing the "Unified Model of Lunar Translation/ 
Rotation" (UMLTR) currently being developed by the author for 
eventual use m the analysis of the Lunar Laser Ranging Experi- 
ment (LURE) data. 

This paper first presents a description of the UMLTR. Secondly, 
numerical results of an investigation of the effects of the 
mutual potential "figure-figure interaction" terms on the physi- 
cal librations are presented. It is shown that these terms 
should be included in a lunar physical libration model when 
accuracies of 0.007 to 0.01 arc seconds (2 cm to 3 cm range pre- 
cision) are required. 


* This work has been supported by the National Aeronautics and 
Soace Administration undeT grant NSG 1152. 

° Present address. 


J D Wulhoiland (ad }. Scientific -ivo hcanons of Lunar Laser Ranging 6i- 7 7 Ad Rights Reserved 
Copyright © 1977 by D Rezdel Publishing Comvany Dordrecht- Holland 
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1. INTRODUCTION 

The ever- mere as zng accuracy of the Lunar Laser Ranging' Experi- 
ment (LURE) data has necessitated a re-examination of the theor- 
etical models used in calculating the lunar orbit and lunar 
rotational motion ♦ 

Previous studies have shown that numerical or M special perturba- 
tion” models axe better suited to provide the required accuracy 
than available analytical models. Thus, numerically integrated 
lunar orbits such as described by Oesterwinter and Cohen (1972) 
and Mulholland and Shelus (1973) have replaced the use of vari- 
ous versions of Brown *s Lunar Theory. Also, solution of the 
lunar rotational equations of motion by numerical integration as 
done by Williams and Slade (1973) and Papo (1973) supplemented 
the use' of analytical theories such as those of Koziel (1948) 
and Eckhardt (1965, 1970, 1973). 

The above numerical models have made the assumptions that (1) 
lunar rotational and translational motions are uncoupled, and 
(2) Earth and lunar radii are small compared to the mutual 
Earth-tMoon distance, r. 

Duboshin (1958) has shown that translational and rotational 
motions of two interacting rigid bodies axe rigorously coupled 
if terms are retained in the differential equations proportional 
to the inverse* fourth power (r~ 4 ) of the mutual distance between 
the bodies. If terms proportional to r~ 4 (r~ 5 ) are retained 
then second (second and third) degree gravitational harmonics 
of both bodies appear m the forces and second and third (second, 
third, and fourth) degree harmonics appear in the torques. Also, 
it can be shown (see section 3) that if terms proportional to 
r~ 5 are retained in the differential equations, then "figure- 
figure interaction” terms from the mutual potential appear in 
the torques exerted on both bodies. 

Existing models, however, indirectly account for coupled lunar 
translational and rotational motions. This is usually done by 
assuming that Cassini's laws are satisfied or by using an analyt- 
ical theory for the lunar rotational motions while calculating 
the orbit. This precomputed lunar orbit is then used as an input 
to the numerically integrated rotational equations. This orbit 
can then be used in a recalculation of the lunar librations, 
etc* 

The second assumption — to the author 1 s knowledge — has not been 
accounted for in existing models. The consequences or assump- 
tion (2) are that the Earth is treated as a particle m the 
field of a non-sohencal Moon and the Moon is treated as a par- 
ticle m the field of a non-spherica 1 Earth m the same model. 
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THE EARTH-MOON MUTUAL POTENTIAL ON LUNAR PHYSICAL LIB RATIO NS 

In other words, the rigorous mutual potential, as discussed by 
Brouwer and Clemence (1961) of the Earth-Moon system is not 
modeled, 

Eckhardt (1973) , JCaula and Baxa (1973), and Pesek (1975) have 
shown that the effects of the Earth acting on the third and 
fourth degree lunar gravity harmonics is observable within the 
current LURE data accuracy (10 cm) . It can also be shown that 
torques acting on the Moon due to the interaction of the Earth's 
oblateness with lunar second degree coefficients are roughly of 
the same magnitude as those produced by the Earth acting on cer- 
tain fourth degree lunar harmonics. 

The preceding considerations, viz,, (1) the trend to numerical 
integration, (2) translation-rotational coupling, and (3) mutual 
potential terms, led to the development of the Unified Model of 
Lunar Translation/Rotation (UMLTR) . This model is essentially 
a numerical integration of the coupled translational-rotational 
motion of the Moon simultaneously with the translational motion 
of the planets. The heliocentric motions of the Moon and planets 
are integrated in a manner similar to that of Oesterwmter and 
Cohen (1972). 

This paper presents a brief description of the UMLTR, It also 
presents a numerical investigation of the effect of previously 
neglected mutual potential terms on the lunar physical librations 
using the UMLTR. 

The rotational motion of the Earth is omitted from these consid- 
erations since the usual practice is to obtain its rotational 
motion from observational data such as the BIH data. This prac- 
tice is to be followed at least m the initial use of the UMLTR. 

2. THE UNIFIED MODEL OF LUNAR TRANSLATION/ROTATION 

The Unified Model of Lunar Translation/Rotation (UMLTR), as 
currently programmed treats tne Sun (1), Mercury (2), Venus (3), 
Mars (6), Jupiter (7), Saturn (8), Uranus (9), Neptune (10), ana 
Pluto (11) as spherical homogeneous particles. It treats the 
Earth (4) and Moon (S) as aspherical rigid bodies. 

The equations of motion for the above system referred to a helio- 
centric frame whose axes (X x 1 } translate with respect to an 
inertial reference frame (mean equator and equinox of 1950.0) 
aTe: 

{i\} * G{(an - 

* {m-> U DP > - g{Z^ "jVlj} (con*.”} 

ih 
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and 


- {m.-l 7 x Uu 3 } T {m-1 ^ - mf* Fl } 

■>* {nu^ 7^ U^ii - nr 1 7^ U^} 

+ Uls + m l 1 ^5 U 15 ^ > (1 = 2 , 11 ) 


CD 

Cconcl 1 d„) 


( V‘T 5 { CtHjHOjH} . 


Cj,k = 0, 1, 2, 3). 
In the above equations: 


(2) 


.DP 


11 11 

= g EE 

l>j=2 


m. m . 
i 3 



the direct planetary terms. 


U 45 


u 14 

«15 


n 

1 

1 


(8 



is the Earth-Moon system mutual potential in the form 
derived by Giacaglia and Jeffreys (1971) , 

is the potential for the interaction of the Sun with 
the asphencal Earth, 

is the potential for the interaction of the Sun with 
-the asphencal Moon, 

includes relativity perturbations on all planets using 
the Eddington/Robertson form as presented by Anderson 
(1974) , and tidal accelerations on Earth and lunar 
orbits formulated similarly to that by Oesterwinter 
and Cohen (1972), 

is an orthogonal matrix of Euler parameters 
CM = 0, 1, 2, 3), and 


<V 


is the angular velocity of the Moon with respect to a 
defined reference frame (j =0, 1, 2, 5). 


The rotational motion of the Earth is adopted from Woolara (1953) 
in the current version of the UMLTR. 


Equations (1) are standard m form except for tne mutual poten- 
tial terms m Uu 3 and the implicit coupling between equations 
(1) and (2). 

Equations (2) are the differential equations for the motion of a 
set of selenocentric — "body fixed"— principal axes, (t }, with 
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THE EARTH-MOON MUTUAL POTENTIAL ON LUNAR PHYSICAL LIB RATIO NS 

respect to a set of reference axes (Z }. The relative orienta- 
tion of these frames is shown m figure 1. 

The orientation of {Z^} with respect to inertial space can be 
completely specified once the relative position of the Earth and 
Moon are known with respect to In figure 1, the axes 

{X x 'T} 

are geocentric axes translating with respect to the frame 
{X x ' }. Thus, r,A,<? are the spherical polar coordinates of the 
Moon relative to the Earth referred to {X^'^}. Finally, the 
unit vectors {Kj.} of the reference^ frame {Zj.} are related to 
the spherical polar unit vectors r,X | via 



Figure 1. Coordinate frames for lunar rotational motion. 
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The relative orientation of Cz-^} and {Z x } is specified by 
the Euler parameters So > 61 , 82 >83 using the notation of Morton 
and Junkins (1973) . Accordingly, 


where 


[c] = 


a.} = [c] {z.} 

0O 2 + Sl 2 -B2 2 -83 2 

2(6iS2-8 0 3 3 ) 

2 (8183+8082) 


C4) 


2 ( 8 i 8 2 + 8 q 8 3 ) 

e 0 2 -Si 2 - 8 2 2 -B 3 2 

2 ( 6283 - 8081 ) 


2 ( 6183 - 8062 ) 

2 ( 8283 - 8081 ) 

8 o 2 ~ 8 i 2 - 6 2 2 + 8 3 2 





of 

equation i 

B o 


-8 2 

-83 

81 

So 

-83 

32 

32 

63 

3 o 

-Si 

S 3 

-s 2 

81 

So 


(5) 


The angular velocity vector {Q.} is essentially the difference 
m absolute angular velocities if the Moon <o 2 m 3 , and the 
reference axis system -X sin d, o, X cos 0 expressed m com- 
ponents in the system {2 }, viz. 


U.} = 


0 


\ 

0 

r 


« 

on 


-X sin <P 

t 

- [C] ■ 

« 

u 2 


9 

f 


* 

^3 


X cos <h 


( 6 ) 


The time derivatives of (Q.) as required m equation (2) can 
be evaluated from Euler's equations for a rigid body, vis. 


* f 
Oil 


' t[/A' ' 

/ i f t 

a 032^3 

• r 

oi 2 

* = < 

t 2 /b» 

- -8'<4»S 

• I 

^3 , 


T 3 /C' J 

[ Y ; 


(7) 


The torque T* (1 - 1, 2, 3) includes terms representing- (1) 
the action of 1 ! spherical Earth and Sun on lunar second degree 
harmonics, (2) the action ox a spherical Earth on lunar third 
and fourth degree harmonics, and (3) the interaction of a tn- 
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axial Earth with lunar second degree harmonics as specified by 
the Earth-Moon mutual potential. 

The set of second order equations (1) and (2) are integrated 
numerically — in the present version of the (JMLTR — using an 
implicit single-sequence method developed by Everhart (1974) , 

The method uses "Gauss -Radau spacings for substeps- within each 
sequence," and integration orders of 7 to 31 are possible. 

The following translational options are currently available m 
the program: (1) integrate all planets, ( 2 ) integrate Mercury, 

Venus, Earth, Moon, Mars, Jupiter, and (3) integrate Earth and 
Moon only. 

In options (2) and (3) above ,the remaining planetary ephemerides 
are read from JPL ephemeris tapes. 

A comparison of the calculated values of the lunar physical 
librations p, la, and r as defined by Eckhardt (1965, 1970), 
with the integrated values produced by the UMLTR is in progress. 
The results of the comparison are to be shown elsewhere, 

Breedlove (1976) . 

3. FOURTH DEGREE TORQUES ON THE MOON RESULTING FROM THE MUTUAL 
POTENTIAL 


The mutual potential of two arbitrarily shaped rigid bodies may 
be put into the following standard form as first shown by 
Giacaglia and Jeffreys (1971) : 

Gmum 5 ° n / I \ n 

u 45 -= -7—- £ E p ™C sin (^) *nm ' C8) 

3 n -2 m =0 ' ' 


where 

$ = X cos ml .45 t Y sin 111 A 45 . 

nm nm nm 

The coordinate system implicit in equation ( 8 ) is centered at 
one of the bodies and r 4 s, 645,145 are th® spherical polar 
coordinates of the mass center of the other body with respect 
to that coordinate system. 


Representative values of the coefficients X nm and Y^ are 
listed here as presented m Giacaglia and Jeffreys (1971) * 


X_ = a 2 C_ + a’ 2 C' , (j = 0,1,2) 

2) 2s 2j 

r 2] " l2S 2j * a ' 2S 2] ’ (J ■ U2) 


(9) 
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X^o " + a^C^Q 

+ Sa^a 1 ^(CzqCiQ * ^ZZ^ZZ + ^22^22^ j 

etc* 

The X^ and Y^ are composed of two - types of terns, viz. 

(1) "particle-figure interaction'* (PFI) terms such as a 2 C 2 j , 
a^C^o, ^a'^C40. and (2) "figure- figure interaction" (FFI) terms 
such as’a?a , 2 C 20 c 20 * equation ( 10 ) a geocentric principal 
axis system is utilized, thus the Cy and Sy are constants. 
Since the Cy and Sy must be evaluated in the same frame, 
they become functions of the relative orientation of the two 
bodies. 

Torque expressions may be derived from equation ( 8 ) by forming 
partial derivatives of U 45 with respect to suitable rotational 
variables such as the direction cosines specifying the relative 
orientation of the two rigid bodies as shown byBeletskii (1966). 

Torques due to "particle- figure interactions" and involving the 
third and fourth degree lunar harmonics may be found in the 
works of Kaula and Baxa (1973), Pesek (1973), and Eckhardt (1973). 

Torques due to "figure- figure interactions" of a triaxial Earth 
with a Moon represented by second degree harmonics have been der- 
ived from equation ( 8 ) and are listed below: 

Tj/A' = Gim l a 2 a ! r _ 5 C2o[ 18 P40T , V' 

+ SPitifO’Y" + Y'a")cA - (T ' 8" - Y , '3 , )sX} 

+ (? 4 2 / 2 ) { (ct 1 a" - S'S")c2A - ( 8 'a" + a’3")s2A}] 

+ GHi li aVr _5 C22[-6P4a(3 , 3" - a '“ ! ') 

+ 3Pui{-(a"Y' + ‘f" a 1 )cA - (Y'8" y"3’)sA> 

+ 3P4 2 {y’y"c2A} 

+ (P 43 / 2 ) { (a"Y’ - 1 - Y"a’)c3A + Cy'3" * Y m S')s5A} 

+ (? 4 .u/ 4 ) {(o' a" - S'3”)c4A 

t (3* a" + a'3")s4A}] , (11) 

T 2 /B 1 = Gmua 2 B , r' 5 C 2 o[-lSPuoYY" 

+ SPuitOY" - yo")cA - C3Y” + Y3”) SA } 
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+ CPu 2 /2MC8B" - aa")c2X - (a?” - ga")s2X}] 

+ Gm ll a 2 g , r" 5 C 2 2[-6P4 0 (S8" - csa") 

+ 3Pt].{(ctT" - ya")cX - (3y” + yS")sX> 

- 3P4 2 (yy"c2X> 

~ (P43/2) { (ay 11 + ya'')c3X + (By 1 ' + yB") s3X) 

* (P44/ 4 HCSB" - aa' f ) c4X - C<s8" * Sa")s4X}] 

( 12 ) 

T’/B' - Gnn 1 a 2 7 , r- S C 20 [18P40YY' 

+ 3P4]_{ (ay ' -- a'y)cX (By 1 A y 6 ')s^} 

+’ (P42/2) { (ao f - 8B')cX + (a8' + 8a')sX}] 

+ Gm4a 2 y'r" 5 C 22 [6P4 0 (aci' - SB') 

+ SPij^-Cay' * ya')cX + (By’ + yS')sX} 

+ 3P42{yy'c2X} 

+ (p43/2){(ay' ya’)c3X + (By' + yB')s3X} 

+- (P44/4) { (act ' - BB')c4X 

(aB‘ + 3a')s4X}] (13) 

where the relative orientation of the geocentnc pnncipal axis 
frame {y^} and the selenocentric principal axis frame (tx) 
is specified by 

a ct' ct" 

{y ± } = S 3’ 3" (2 i > (14) 

_ Y Y 1 Y" 

using the notation similar to Beletskii (1966), In equations 
(11) to (13) the subscripts <15 have been omitted for clarity 
and the abbreviations c = cos and s = sin have been used. 

The torque expressions in equations (11) to (13) contain the 
direction cosines a,y',y" of the rotation axis of the Earth 
with respect to the selenocentric frame (ti). It can be shown 
that these quantities exhibit periodic variations with the lunar 
nodal period (18,6 years), among others. 
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Lure (1963) presented a vector-dyadic form of equations (11) to 
(13) when the interaction of the Earth's oblateness with a tri- 
axial Moon was considered. 


4. THE EFFECT OF FIGURE-FIGURE INTERACTION TERMS OF THE EARTH- 
MOON MUTUAL POTENTIAL ON LUNAR PHYSICAL LIBRATIONS 

A qualitative analysis of the effect of the "figure- figure inter- 
action" (FFI) torques on the lunar physical librations was under- 
taken in the following manner. The UMLTR was used to produce a 
simulated set of physical librations over a three-year period 
beginning at J.D. 2441200.5. This set ("observations") was 
calculated including the effect of the usual second degree 
terms m the lunar torques as well as the FFI torques due to the 
Earth's oblateness. A "calculated" set of librations consti- 
tuting the model was next produced in a similar manner but with 
the FFI torques omitted. The model was then fitted to the obser- 
vations in a least squares sense by adjusting the lunar inertia 
ratios a ' , 3 1 , and y 1 and the rotational initial conditions 
( 6^)0 and (Si)o- The results of this fit dre shown in 
figure 2 . 

The nature of the residuals indicates that the FFI torques pro- 
duce monthly variations m the latitude librations, 0 and la, 
with amplitudes of about' .002". Also, the residuals m p, la, 
and t indicate a long period effect that could not be fit over 
the adopted time interval (3 years) . This long period effect is 
due to the 18.6 year variation m the direction cosines y,y' ,y" 
of the Earth's spin axis with respect to the lunar frame (z } 
as discussed in section 3. Physical librations m p, la, and 
t with amplitudes of 0 . 01 " are apparent from figure 2 due to 
the presence of the FFI torques. 

The above results were obtained with nominal values of 
6 = 0.6293 x 10 “ 3 , y = 0.227 * 10 “ 3 , and‘C 2 o = -1.08261 * 10 ~ 3 . 
Data was obtained in 5-day increments and the 11th order version 
of Everhart's (1974) integrator was used. Planetary and lunar 
initial conditions were adopted from Oesterwmter and Cohen (1972) 

5. CONCLUSIONS 

The Unified Model of Lunar Translation and Rotation has been 
developed for studying lunar translational/rotational coupling 
and the effect of terms arising m the Earth-Moon mutual poten- 
tial. It is being verified by comparisons with Eckhardt (1965, 
1970, 1973, 1976) and Oesterwmter" and Cohen (1972) for eventual 
use in the analysis of LURE data. 

Figure-Figure Interaction torques produce effects on the lunar 
physical librations comparable in magnitude to those produced by 
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certain fourth degree terms m the lunar gravity field. In par- 
ticular, monthly variations of amplitude .002" were found m 
the latitude lib rations, p and la. Also, long period variations 
in p and la and the longitude librations, t, with amplitudes 
of -0,01" and periods that are fractions and multiples of the 
lunar nodal- period (18.6 years) were using a numerical experi- 
ment, The larger of the above effects causes range variations 
of "10 cm. Thus, the FEI torques should be included in a lunar 
libration model accurate to "0,01". 


6 . SYMBOLS 


ORIGINAL 
0 F POOR QUAU'T*’ 


All primed quantities refer to the Moon, unprimed quantities 
refer to the Earth. The following symbols have not been 
defined previously: 


r. 

i 

m. 

i 

a’ 

8' 

r' 

a 



t 

P 

G 


heliocentric position vector to planet 1 
mass of planet i 
= (C f - B^/A 1 

- (C 1 - AO/B 1 

- (B 1 - A 1 ) /C 1 
reference radius 

gravity harmonic coefficients 
gravity harmonic coefficients 
physical librations in latitude 
physical libration m longitude 
Associated Legendre polynomials 
universal gravitational constant 
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DISCUSSION 

Counsetman; Have you compared your theory with lunar laser 
observations yet? 

Breedlove: No, we have not reached that stage, although we 

have the data. 

Mulholland: I think that the reason for your large residuals in 

the comparison of Vibration models is probably due to having 
chosen a rather poor analytical model against which to com- 
pare. The errors in the comparison staandard will give spurious 
corrections to the numerical starting conditions. Since the 
computer "Moon" has an infinite Q, this will stimulate 
artificial free Vibrations in the' numerical integration, You 
must use a better standard. 

Breedlove: I have recently received the improved Eckhardt 

series 300, which I will use m the future. 

King: Eckhardt now has an even newer theory, designated 
series 400. 
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A Unified Model of Lunar Translation/Rotation 


Dr. William J. Breedlove, Jr. 

Associate Professor 

Mechanical Engineering and Mechanics Department 
Old Dominion University, Norfolk, Va. 23508, USA 


ABSTRACT: 

A unified model of lunar translation/rotation (UMLTR) is 
described in which the heliocentric motions of the planets, the 
geocentric motion of the Moon, and the rotational motion of the 
Moon are integrated simultaneously. Novel features of the UMLTR 
include l) inclusion of translational/rotational coupling, 2) inclu- 
sion of figure-figure interactions (FFI) between the Earth and Moon, 
3) the use of Euler parameters for Lunar orientation, and 4) the 
use of a relatively new integrator developed by E. Everhart. At 
present the UMLTR fits Eckhardt’s 300-301 series with standard 
deviations of 0.1" , 0.15", and 0.1" in p, 1^, and x. The UMLTR 
• lunar orbit differs from that in DE96 by 80 meters in 3 years and 
200 meters in 6 years. Coupling effect residuals of 5 cm. ,,0 007" , 
and 0.001" after 3 years arise in lunar semi -major axis, perigee, 
and longitude respectively if fourth order lunar harmonics are 
included over and above a third order model. 
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1. INTRODUCTION AND SUMMARY 


The Lunar Laser Ranging Experiment (LURE) has placed three lunar 

ranging retro-reflectors on the Moon as described by Bender (1973). 

A series of measurements of the distance to these retro-reflectors 

has been inferred from measured round trip travel times of pulsed 

laser beams. These measurements have been made regularly since 

August 1969 and allow determination of the distance to about +_ 8 cm. 

as reported by Williams ( 1976 b) . This same report indicates that a 

resolution of +_ 2-3 cm. is expected eventually. 

> 

According to the LURE team, the LURE data can be used to obtain 
improved knowledge of 1) retro-reflector and earth-based telescope 
coordinates, 2) lunar center of mass motion, 3) lunar librational 
motion, 4) Earth rotational motion, 5) Earth continental drift, 

6 ) gravitational theories, and 7 ) the internal composition of the 
Earth and Moon. Indeed, a number of major announcements have 
recently been made by the LURE team including l) a confirmation of 
the equivalence principle of general relativity reported by 
Williams ( 1976 a), and 2) confirmation that lunar laser ranging is 
capable of accurately determining the Earth’s rotation as reported 
by Stolz (1976). 

The ultimate accuracy of the LURE data has led a number of 
investigators to look into previously neglected physical and 
dynamical effects in existing models for the translational/ 
rotational motions of the Moon. Thus Gurevich (1976) analyzed the 
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tidal effect in the theory of lunar rotation; Calarne (1976) has 

I 

studied the presence of small free lihrations; Lambeck (1975) 
has investigated the effect of tidal dissipation in the oceans 
on the Moon’s orbit; and Melosh (1975) has studied the effect of 
Mascons on the Moon's orientation. 

Existing models for the translational and rotational motion 
of the Moon generally contain the assumptions that (l) lunar 
rotational and translational motion are uncoupled, and (2) Earth 
and lunar radii are- small compared to the Earth-Moon distance. 

A special, perturbation model not incorporating the above tvo 
assumptions is currently being developed by the author. This 
model, the Unified Model of Lunar Translation/Rotation (UMLTR), 
is being used to investigate the above assumptions. This' paper 
documents the UMLTR, provides comparisons with existing models, and 
gives the results of certain numerical experiments. 

Preliminary results using an early version of the UMLTR were 
presented by Breedlove (1976a) showing that figure-figure inter- 
actions (FFI) of the Earth and Moon have at least a 2-3 cm. effect 
on the distance measured to the retro-reflector. These PPI can 
be investigated if the above assumptions are not made. 

Novel features of the UMLTR include l) the ability to integrate 
the coupled rotational -translational motion of the Moon, 2) inclusion 
of FFI torques based on the Earth-Moon mutual potential, 3) the lose 
of Euler parameters for lunar orientation, *0 the use of a relatively 
new integrator developed by Everhart (197*0; and 5) the option to 
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integrate all planets, or only the inner planets or only the Earth- 
Moon barycenter together with the lunar translational -rotational 
motion. In the latter two options, DE 96 provides the remaining 
planetary positions and velocities. 

Section 2 provides a summary of the equations of motion with 
a brief description of the numerical integration scheme used. 

Section 3 provides the results of certain verification tests of 
the UMLTR including a comparison with DE 96 for 6 years for transla- 
tional motion; a closure test for 6,5 years of integration as a check 
on numerical roundoff errors and a 3 year comparison of the UMLTR 
with Eckhardt's (1976) 300-301 series for the rotational motion. 

Section k provides the results of some numerical experiments involving 
coupling and the effects of FPI. 

At present, the UMLTR fits Eckhardts 300-301 series with standard 
deviations of 0.1”, 0.15", 0.1" in p, 1^, and t. The UMLTR lunar 
orbit and that of DE96 differ by about 80 m after 3 years of integration 
and 200 m after 6 years of integration. 

The effect of coupling is shown indirectly by plotting the lunar 
orbital residuals between UMLTR runs which differed only by the 
assumed lunar mass distribution. After 3 years, differences in semi- 
major axis, perigee, and longitude are 5 cm. , ,0 007", and 0.001” 
respectively when a run containing Uth order lunar harmonics is 
differenced with one containing 3rd order harmonics. 
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2. EQUATIONS OF MOTION 


2.1 General Equations 

The equations of motion for the Moon and planets used in xhe 
present version of the UMLTR are summarized here. The heliocentric 

A 

motions of the planets and the Earth-Moon barycenter are integrated 
simultaneously with the geocentric translational motion of the Moon 
and its rotational motion. 

The equations actually programmed are based on the more general 
equations: 


' m.r . = V.U + F. , (l) 

{ V = !af<V { V } > (2) 

with i = 1 , . . . ,11 ; j , k - 0 ,1 ,2 ,3 . 

In equations (l) and (2), r. is the position vector of the ith 

— i 

body with respect to the solar system barycenter and l J are 

Euler parameters providing the orientation of 'the Moon’s principal 
axes with respect to a reference axis system defined later 
(section 2 . 6 ), 

The potential function U in equation (l) is of the form 

11 11 11 , 

U = G 2 Z m.m,rT7 +' Gnu Z m.r7. 
i-2 J=2 1 3 13 1 j=2 j 

irj 




© 
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( 3 ) 


+ + 2 U lp 

p=2 

= u + + E u lp 

p=2 


where the function U ap appearing in this equation represent the- 
mutual potential between bodies a and p given in the form 
established by Giacaglia and Jeffreys (1971 ) viz-. 


u ap = 


Gm m 

a p 


CO n 
2 2 
n=2 m=0 


p (s in 


nm 


ap' crp ap 


(M 


where 


_ = X<JP 
ap run 


cos m 



- Y 00 

nm 


sin m X_ . 

ap 


The quantities X° P and are functions of the harmonic 
nm • nm 

coefficients of bodies 0 and p and r , 4 , X are the sphencal 

K ap’ Y ap’ ap * 

polar coordinates of one body with respect to a coordinate system 
at the first body- The notation used for the "N" bodies is Sun = 1, 
Mercury = 2, Venus = 3, Earth = U , Moon = 5» Mars = 6, Jupiter = 7> 
Saturn = 8, Uranus = 9, Neptune =10, and Pluto = 11 . 

The force vector F. is given by 


F. 



Ff + F t 


(5) 


2* 

where 5\ represents the relativistic perturbations based on an 
isotropic, Parameterized Post-Newtonian N-body metric and the 


27 



vector F\ represents the perturbations produced by the asteroids. 

I* 3# 

Expressions for F. and F. vere taken from the discussion of the JPL 

Development Ephemeris Number 96 by Standish et. al. (1976). The 

vector represents the tidal coupling perturbations and was 

modeled as suggested by Oesterwinter and Cohen (1972). The vector, 

{$.}, and matrix, [B ] in equations (2) are composed of the Euler 
i kj 

parameters 8 Q , B 2 , 

The angular velocity vector {Q.} in equations (2) is essentially 
the difference between the absolute angular velocity components of 
the moon and a reference axis system as described later in section 
2 . 6 . 


2.2 UMLTR Model for Translational Motion of Planets 

In the UMLTR, the heliocentric motions of the planets except 
the Earth are integrated in the form 

£ 1 = r 1 -? 1 (i = 2 11 , i * 1,5) «5) 


where 


r x = m^V^U + U 1 ^ + U 15 ) 


and 
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mT^V. (U) + m 7 1 F r + mT 1 ? 3- 
l—i l—a 1 — 1 ' 

—1 r ^ 

In equations (6), a simplified version of nu |\ is used. The 

version programmed is the Eddmgton-Eobertson form as given by 

Anderson ( 1974 ), viz.. 



m.F^f 

i— i 


G m 1 c** 2 p i 3 { [2(0+7)^ - yV 2 ]^ 


+ 2 (l+Y)(P i P i )P i ) 


(7) 


•where 


^ = % *£i 
♦i = 


and 0 and y are relativity parameters. Note that the term 
is not included in equations (6) which, according to Moyer (1968) 
is a reasonable omission for accuracies of 1 m for the inner planets 
and 1 km for the outer planets. 

The effect of the asteroids -was included only as they perturbed 
Mars orbit. Standard orbits for the asteroids was adopted as given 
in Standish ( 1976 ). 

In equations (6), the Sun is treated as suherical; and C„ 

20 30 

1 1 i 

are retained for the Earth; and C 2Q , C 22 , and S 22 are retained for 
the Moon. 
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2.3 UMLTR Model for Motion of Earth-Moon Barycenter 


The heliocentric motion of the Earth-Moon barycenter in 
UMLTR is integrated in the form 

% = m 5 (m u + n^r 1 £_ + + n^)" 1 

where 

4 = Ik - Si' 

Using equations (l ) , can be written as 
% = -G(m 1 +m u +m 5 )[m 4 M" 1 £ tt p'J + 

+ M ’ 1 ,L w 3 Cm 5 (e d-% )p 5j " 

0r^»5 

" Mp^p‘3] + + F^] - Mm^ 

+ (m^ + m^ + V^U^^ + ^u 14 ] 


where p. = Gm 

O J 


and 


M = m^ + m 5 . 


In equation (9), m^ ME^ is omitted and the remaining 
relativity terms come from equations (7) with 


h = Gm l p lJ + Gm 5 p U5 


the 


( 8 ) 


(9) 


.(10a) 
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and 


<j> 5 = Gn^ + Gm^p“^. (10b) 

2.4 UMLTR Model for Geocentric Lunar Motion 

The geocentric lunar motion of the Moon in the UMLTR is 
integrated in the form 


*« i« 


A} % ~ £4 ’ 


( 11 ) 


Again using equations (l), can written as 


% = " G £ M P’ 3 " g “ 1 [%P5° " 


-3 


-3- 


11 


-3 


-0 £ .jKfij - £j)p 5j 

0^,5 


+ [VjD 15 + 


(e^-Rj^uP 


- % lSi,v Lh + £ + eJi 


£4 T £4^ 


+ Mm ^ 1 ^ 5 - v/ 5 ], 


( 12 ) 


In equations (12), U 4 ^ retains C^q and C^q for the Earth and 

1 t 1 

Cgo* G 22* ^22 ** or ^ oon * < ^ 1S coupling terms are 

computed as suggested by Oestervinter and Cohen (1972) viz.. 
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■V £ - (M,. *'j 'Ti? - 1 ' 5 “i4 


(fir x fip) x fir 

— -- - (13) 


l k ' m 5' 


1% x 


vhere 6. . is the delta function and C is the tidal coupling 
coefficient. The relativistic terns in equations (12) are computed 
as suggested hy Moyer (1968) as 

” m b^. = ^M (S) " 4 (S) 


+ %(=) - 4 (M) 


(lb) 


vhere r,.(S) means the relativistic acceleration of the Moon due 
“M 

to the Sun, etc. In particular r^(S) and (S ) are calculated 

from equations (7) and (10). The terms r.,(E) and r_,(M) are 

~Vi 

calculated using the "N-hody" metric in the form given "by 
Standish (1976) viz., 


= 2 v (r - r. )r7?{- 2 (B+y)c~ 2 £ y r^ 

1 j#i J J 10 k#i 151 “ 


{ 2(3-1) c" 2 Z P.r J + Y(V./c f + (l+Y)(V./c) 2 
^ k ok i 0 


- 2(1+y)c“ 2 ‘ r . - (3/2)c" 2 [(r i -r j ) • 


- 2 , 


- 1.2 


+ (l/2)c~ 2 (r.-r. ) • £.} + c“ 2 £ y,r~^{(r^-r ] 




• • 


• [(2 + 2y)r. - (l + 2Y)r.]}(r. - r.) 

“i “0 “ 1 — 3 
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+ (3/2 + 2y)c~ 2 I U,£.r"^ (15) 

j J 1J 

where instead of summing over " j" , M j M is set equal to 4 for r^(E) 
and "j" is set equal to 5 for r_(M). 

2.5 UMLTR Model for Lunar Rotational Motion 

The rotational motion of a rigid Moon satisfies Liouville's 
equations. A special fom of these equations are referred to as 
Euler's equations, viz., 




which represent a rigid rotation from a set of orhitally defined 

reference axes {Z.} to the lunar principal axes {z.}. These 
1 1 

axes are described in section 2.6. 


33 



The relative orientation of {Z^} and {z^} is given by 

{z.} = [C]{Z.} (18) 

where 

S 0 + B r B 2 " B 3 2 ^ B 1 B 2 + B 0 B 3^ 2 ^ B X B 3 ~ B 0 B 2^ 

[C] = 2(6^2 - B 0 B 3 ) b* - el * e| - Bj 2(B 2 B 3 + B^) 

_ 2 (B 1 B 3 + B 0 B 2 ) 2(B 2 B 3 - B 0 8 1 ) B* - &l - B 2 + B§_ 

using the notation of Morton, et al. (197^ )♦ 

The absolute angular velocity of the Moon, [ ] is given by 
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In equations (19) > [C] is matrix [C] augmented, by a row and 

column of zeroes - the zeroes forming the first row and column 

of [c]. Also, in these equations the angular velocity of the 

* * * 

reference frame, {z^}, is given by l-Asin<j>,<j>,Acos4>j . Finally if 

equations (19) are solved for {$.} and put in a second order 

J 

form, then the lunar rotational equations assume the form given 
in equations (2) viz. , 


where 




<v 


r -\ 


r 


0 


0 

w i 


-Asin<j> 


> - t°i < 

. ) 

W 2 

$ 

W 3 


♦ 

Acos<{> 

J 



J 


The time derivatives of {ft.} as required in equations (2) are 

J 

evaluated from Euler's equations, (l6). 


2.6 Coordinate Frames 

The basic inertial reference frame, X_^ (i=l,2,3), is located at 
the Earth's mass center and is the traditional mean equator and 
equinox of 1950.0 system. Axes x^ , and y^, with unit vectors and 
are Earth and Moon principal axis systems. Axes Z^, with unit 
vectors K. , are the reference axes for the lunar rotational motion. 
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The unit vectors K. have a defined relation to the unit vectors 

- i 

ij .5 a spherical polar system locating the Moon with 

respect to , viz . , 

-1 = “ir* -2 = "-A* -3 = 

2.7 UMLTR Model for Earth Rotational Equations of Motion 
Although the rotational equations of motion for the Earth 

have "been derived and coded in a manner similar to that done for 
the Moon (Section 2.4 ) those equations are not integrated in the 
current version of the UMLTR. Instead the nutation, precession, 
and spin of the Earth are defined "by the formulae given hy 
Lundquist and Veis ( 1966 ). 

It was found necessary to use the rigorous rotation matrix 
for nutation however, rather than the linearized version incorporated 
there. The UMLTR can also read the JPL DE 96 Ephemeris tape for 
nutation and nutation rates. 

2.8 Forces and Torques 

The forces utilized in the UMLTR have been outlined in preceding 
sections. Detailed derivations are presented in Breedlove (1975 > 
1976h, 1977 ) • The , lunar torques coded in the UMLTR are listed in 
Table 1. Pesek (1973) did not list all of the terms in the 4th 
order particle- figure interaction (PFI) torques and made some 
approximations in the remaining terms. The full fourth order 
PFI torques were coded in the UMLTR, however. A convenient vector 
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form for the FFI torques when only C^g is considered is given by 
Chohotov (1964) in the form 

T = -(15/2)u 1i 0 20 fV 5 [(1 - 7 cos 2 e)i r • I X 1^ (20) 

+ 2 cos 9(i 3 • ! x ^ ^ • I x i 3 ) - | ^ • I X ls ) 


where 8 = * i^. is i' unar ineruia dyadic. 


LUNAR 

GRAVITY 

COEFFICIENTS 

INVOLVED 

TYPE 

OF 

INTER- 

ACTION* 

BODIES 

INVOLVED 

REFERENCE 

2nd Order 

PFI 

Earth , Sun 
Moon 

Eckhardt (1970) 

3rd Order 

PFI 

Earth, Moon 

Pesek (1973) 

4th Order 

PFI 

Earch, Moon 

Pe^ek (1973) 

4th Order 

FFI 

Earth (C , 

C ) dU 

^22 , 

Moon (C 20 , 
°22 J 

Breedlove- (1976a) 

i 

*PFI = Particle - Figure Interaction 
*FFI = Figure - Figure Interaction 



TABLE 1. LUNAR TORQUES 


2.9 Numerical Integration Method 

A relatively new integrator, as described by Everhart (1974a, 
1974b) , is being utilized m the UMLTR. This integrator uses an 
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implicit single sequence algorithm. Substeps within each sequence 
are taken using Gauss-RacLau spacings for increased efficiency and 
accuracy. The integrator has the form of an empirical time series 
method and has been shown to be related to the implicit Runge— 

Kutta method. Orders of 7 >11 ,15 ,19 ,23 ,2? and 31 may be chosen and 
the 11th order option is being utilized on all UMLTR runs. 

3. NUMERICAL VERIFICATION TESTS 

Results from numerical experiments and from several comparisons 
of the UMLTR with existing numerical/analytical models are presented 
here. All runs described here were made on the Old Dominion 
University DEC-10 computer in double precision. This computer 

f 

maintains a precision of 18 decimal digits in double precision with 
72 bits per double precision word. 

3.1 Comparison with DE 96 

Various comparison runs have been made with DE 96 for the purpose 

of validating the UMLTR for integration of the Moon and the planets. 

in the run described here, planetary initial conditions were 

obtained from DE 96 for Julian Date 2^1200.5 and all constants used 

in the DE 96 run were duplicated in the UMLTR. The initial conditions 
* * 

for fL, (i = 0,1, 2, 3) were obtained from Eckhardt’s (1976) 300 
series. All rotational, physical parameters were set to those values 
specified by the 300 series. A program, EMAN, was developed to 
(in part) convert the physical librations of the 300 series into 
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the Euler parameters 8^. Initial rates 8^ are calculated by 
numerical differentiation of a cubic spline fit to 8^ in EMM. 

The basic differences between DE96 and the UMLTR in terms 
of the planetary integrations are that 1) the full relativistic 
equations as given by Standish (1976 ) have not been programed 
in the UMLTR, (see Section 2.2) and 2) the UMLTR integrations 
assume the Sun's mass center is the origin of an inertial 
reference frame. With the above initial conditions, a run was made 
for six years. The results of the comparison between the UMLTR 
and DE96 are presented in Table 2. The quantity tabulated 
there is 


5 r i^DE96 X UMLTR^ + ^ 


DE96 


_ y )^ + 

UMLTR ' 


^296 ~ \lMLTR^ 


in meters. The differences, 6^, are periodic with Jupiter having 
the largest values probably due to the omission of the Sun's 
acceleration in the UMLTR. An earlier comparison run was made 
without modeling the asteroids Juno, Pallas, and Vesta resulting 
in <$ r 's of 398 m and ^608 m for Mars after IjOO and 800 days 
respectively. Inclusion of the effect of these asteroids on Mars 
orbit in the UMLTR produces the d^'s quoted in Table 2. The run 
described here used the 11th order integrator with a constant 
sequence size of 1 day. CPU' time on the DEC10 for the six year 
run was approximately 102 minutes. 
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PLANET 

Dffi PAST EPOCH 

1*00 d 

800 d 

1200 d 

1600 d 

2000 d 

2400 d 

MERCURY 

56 

115 

54 

179 

190 

245 

VENUS 

102 

377 

l84 

342 

263 

669 ' 

EARTH-MOON 

BARYCENTER 

306 

U 58 

1*93 

391 

152 

196 

MOON (GEOCENTRIC) 

37 

82 

35 

120 

199 

88 

MARS 

162 

873 

82 

641 

752 

268 

JUPITER 

6l 

327 

1,278 

5,202 

7.847 

12,400 

SATURN 

mm 

161 

338 

579 

923 

1,421 

URANUS 

9 

38 

141* 

137 

229 

329 

NEPTUNE 

3 

11 

26 

35 

63 

59 

PLUTO 

150 

13 

13 

50 

79 

144 


TABLE 2. COMPARISON OF UMLTR AND DE 96 FOR SIX YEARS (5 in meters) 



































































The primary goal of this work however , is to produce an 
accurate integration of lunar translational/rotational motion. 

The accuracy of comparison of the Lunar orbit with DE9 6 is 
encouraging although no detailed comparison of the UMLTR and the 
lunar ephemeris in DE 96 (LE44) has been attempted. Also it 
should be noted that options exist in the UMLTR to integrate 
only the Earth and Moon or only the inner planets and Moon while 
reading all remaining planetary positions from DE 96 . 

3.2 Closure Test 

A closure test was performed using the UMLTR as a rough check 
on roundoff error. A run was initiated at Julian Date (JD) 
2441200.5 and terminated at JD 2442400.5. The terminal conditions 
were then used to perform a backward integration terminating at 
JD = 2441200.5. The same initial conditions used in the previous 
test were used here. The RMS position vector error for the Moon 
was 18 cm. and the RMS velocity vector error was 7 cm. /day. In 
terms of the physical librations, the error in p was 0.004" , 
the error in I was 0.07", and the error in T was 0.07" • These 
are an indication of the roundoff errors after about 6. 5 years of 
integration. 

3.3 Comparison with Eckhardt ' s 300 Series 

An iterative weighted least squares technique was used to 
fit the lunar physical librations p, I » and T calculated by the 
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UMLTR with those available from Eekhardts (1976) 300 series. 
Eckhardt's 300 solution was supplemented by a set of '’additive" 
and "planetary" terms supplied by Williams (1976c) and by the 
second order term described by Williams, et. al (1973). Also 
Eckhardt's 301 series terms were utilized to provide the effects 
of Hh degree terms in the lunar gravitational potential. 

Differences of the UMLTR and Eckhardt's solution for p, 1^, 
and T were calculated at 3 day increments for 1200 days. The 
quantity 


2 2 
? ^ P 300 ~ PlJMLTR^i + ^ ICf 300 “ Ia UMLTR^i 

+ (T - T ) 2 

^ 300 UMLTR ; i 

was chosen as the quantity to be minimized where the sum is over 

all "observations". The fit was made by adjusting the initial 

conditions, (0.) and (0.) (i=0,l,2,3) and 9 bias parameters 

1 0 1 0 

using numerical partial derivatives. 

♦ 

The quantities (0.) and (0.) are not all independent so 
1 0 1 0 

that in performing the estimation and in calculating the partial 
derivatives the following constraints are maintained 

+ P 1 + 6 2 + 3 3 = ( 20 ) 
^ 0^0 + Vl + ^ 2^2 + ® 3^3 = °* 
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In practice only 0g 5 0^ > 0g , 3 q 0, 0 2 were estimated with 0^ an ^ 
0 ^ "being calculated from the constraints by the method given in 
Moyer (1971). 

In calculating the numerical partial derivatives with respect 

to (0.) , the usual normalization procedure used when integrating 
1 0 

Euler parameters must be suspended. Usually, these parameters 
are normalized by a factor, QN, where 

<*-$ + <%+<$ + 6 3 

during the numerical integration. 

The bias parameters, b^ (i = 1 ,..., 6 ) are of the form 

P UMLTR = P CALCULATED + \ + V 1 
Ia UMLTR = Ia UMLTR + ^3 + V* 

X UMLTR = T UMLTR + ^5 + 1 

while the remaining bias parameters b^ , bg, b^ represent arbitrary 
constant rigid body Euler rotations about axes {z^}. 

The epoch for the fit was 24^1200. 5 and all initial conditions 
and parameters are consistent with DE 96 and Eckhardt's 300 series 
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for that date. Again initial rotational conditions were calculated 
by program EMAN. 

The results of xhis fit are shown in Figure 1, where residuals 
(0-C) in p, la, and x are shown in arc seconds plotted against 
days past 2441200.5- The means and standard deviations of these 
data are given in Table 3. The fit changed the nominal initial 
conditions at most in the 13th decimal place. The estimated values 
of the bias parameters are given in Table 4. All bias values were 
given nominal values of zero. 


LIBRATION ‘ 

MEAN 

(Arc sec) 

STANDARD 
DEVIATION 
(Arc sec) 

P 

-0.1011D-4 

0.093 

la 

0.1850D-4 

0.148 

T 

0.3602D-3 

0.053 


TABLE 3. STATISTICS OF UMLTR FIT TO ECKEARDT 300- 

301 SERIES 
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Figure 1.- UMLTR fit to Eckhardt 300-301 series. 



4 ^ 



Parameter 

Estimated Value 


0.199” 

"a 

-O.OOV'/d 

b 3 

6.042" 

\ 

-0.011'Vd 

"5 

-0.035” 

*6 

0.292X10"^ ”/d 

*7 

0.010” 

OO 

& 

0.013” 

b_ 

0 . 018 ” 

9 



TABLE k. ESTIMATED BIAS PARAMETERS FOR UMLTR FIT TO 

300-301 SERIES. 

The residual patterns of Figure 1 shov the existence of long 
period free lihrations in p and Id that could not he fit; over this 
data arc. The residuals in T shov that the solutions differ hy 
annual terms of about 0.1 arc second and that the free libration 
in T has been suppressed. The estimated values of the bias 
parameters b^-b^ indicates that the coordinate reference frames 
of the 300 series and the UMLTR don't coincide precisely. 

4. TRANSLATIONAL /ROTATIONAL COUPLING AND FIGURE-FIGURE INTERACTIONS 

Several runs of the UMLTR have been made to provide an indication 
of the effect of the earth-Moon mutual potential on lunar physical 
lihrations and the effect of considering translational /rotational 
coupling for the Moon, 
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4.1 Coupling of Lunar Orbital /Rotational Motions 

Duboshin (1958) has shown that if terms factored by are 

retained in the differential equations of motion that coupling is 

present and that the force expressions contain terms involving 

second degree gravity harmonics while the torque expressions contain 

terms involving second and third degree gravity harmonics. The 

_k _c 

UMLTR models terms factored by in the forces and by p^ in 
the torques. 

A nominal run of the UMLTR was made modeling second and third 
order lunar gravity harmonics. Translational initial conditions came 
from DE 96 at JD = 2441200.5* Rotational initial conditions 
vere determined from a fit of the UMLTR to Eckhardt ' s 300 series 
to eliminate large free librations in T. 

Two additional runs were made using the same translational 
initial conditions but changing the lunar gravity model and the 
rotational initial conditions. In the first case lunar second, 
third, and fourth degree harmonics were retained while in the 
second case only second degree harmonics were retained. 

Rotational initial conditions , in each case , were found from fits 
of the UMLTR to Eckhardt' s series. 

Figure 2 shows the residuals between the nominal and the first 
case run. Residuals in semi-major axis in meters, eccentricity, 
inclination, argument of perigee, node, and true longitude are 
shown. All angular residuals are in seconds of arc. The elements 
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Figure 2.- (Cont.) Lunar orbital residuals between 4th and 3rd degree 

torque models . 
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are vith respect to the fixed 1950,0 system. The major effects 
are the growth in longitude amounting to 0.001” in 3 years, the 
periodic grovth in perigee amounting to ,0 GOT" in 3 years , and 
the grovth in semimajor axis of 5 cm. in 3 years. Identical 
residual patterns were obtained for the case two run but the magni- 
tudes were much greater. After 3 years, the residuals in senimajor 
axis, perigee, and longitude for the latter case are 6.6 meters, 
.09”, and .IT” respectively. Apparent secular motions may also be 
noted in inclination, node, and longitude. The ratio of the effects 
noted above between cases one and two are the same as the ratio 
of the constant longitude bias between the two cases. 

The model in the UMLTR directly includes the effects of 2nd 
order lunar harmonics in the forces governing translational motion. 
The inclusion of 3rd and 4th degree terms in the torques can only 
affect the translational motion through the coupling mechanism. 

That is, the presence of these torques affects the lunar orientation 
which, in turn, affects the forces. 

-k~,2 The Effect of FFI on Lunar Physical Librations 

The effect of FFI -interactions on physical librations was shown 
in a preliminary study in Breedlove (l9T6a). There it was shown 
that these terms have at least an 0.01" effect on lunar physical 
librations and thus should be included m any model analysing the 
LURE data. Yoder (19TT) has subsequently derived approximate 
analytical expressions for the FFI effects on lunar physical 
librations. That study indicates "that the largest terms produce an 
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18.6 year llbration of magnitude .007" in longitude and 0.03" 
in latitude as seen from an inertial frame." 
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NOTATION 


Seconds of arc 


m D 

c 

C 

ct,B,y 
c ij ' s ij 

t i 

c „ ,s. . 
1J II 


Mass of jth particle 

Speed of Light 

Tidal coupling coefficient 

Lunar Principle Inertia Ratios 

Earth Gravity Harmonic Coefficients 

Moon Gravity Harmonic Coefficients 
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ON LUNAR ORBITAL-ROTATIONAL 
COUPLING AND FIGURE-FIGURE 
lilTE RACTIONS IN THE EARTH-MOON 
SYSTEM 

W. J. Breedlove, Jr f * 

Department of Mechanical Engineering and Mechanics 
Ola Dominion University 
Norfolk, Virginia 


Abstract 

A numerical model - the Unified Model of 
Lunar Translation/Rotation (UMLTR) - nas neen 
developed by the author. The model numerically 
integrates the heliocentric motion of ail planets, 
the Earth-Moon baryc enter , and the geocentric . 
motion of the Moon coupled with its rotaticnal 
motion. Results are presented of numerical exper- 
iments performed with the UMLTR whicn show 
coupling - effect ormtal residuals of 5cm. , 
,0007”, and 0.001” m semi-ma^or axis, perigee 
and longitude - after three years - between a 
third-order lunar gravity model and a fourth- 
order model. Figure- figure interactions provide 
perturbations of at least 0,01" in the lunar 
pnysical iibrations. Finally, the possibility 
of estimating simultaneously lunar inertia 
ratios, a, $, y with lunar coefficients c ^ and 

C 22 is discussed based on the UMLTR and the Lunar 
Laser Ranging Experiment (LURE) observable. 


Nomenclature 


A,B,C 

C 2Q ,C 22' S 22 

C 20' C 22' C 30 

M 

R 


Lunar principal moments of inertia 

Lunar Harmonic Gravity Coefficients 
Earth Harmonic Gravity Coefficients 

Lunar mass 
Mean lunar radius 


1. Introduction and Summary 

Consideration of the coupling between lunar 
orbital and rotational motions and tne effect of 
the mutual potential oe tween the Earth and Moon 
have been - for the most part - ignored in many 
numerical and analytical models developed m tne 
past. The complexity of any model incorporating 
these effects and the level of observational 
accuracy available in tne past did not warrant 
their treatment. Recently, however, the observa- 
tional accuracy has improved to such an extent 
tnat inclusion of the aoove effects may be 
important. The Lunar Laser Ranging Experiment 
(LURE) is currently providing measurements of 
the equivalent distance from an Earth-bound tele- 
scope to the lunar retro-reflectors accurate to 
13cm (over the period 1969-1977) . The antici- 
pated accuracy level for tms data is 2-3 cm. 


The author has developed a numerical model 
for the coupled translational /rotational motion 
of the Moon referred to as the Unified Model of 
Lunar Translation/Rotation (UMLTR) , This model 
consists of tne direct numerical integration of 
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and Mecnanics 
Member AIAA 


the heliocentric orbits of all planets, the Sarch- 
Mcon baryc enter, and the geocentric orbit of tne 
Moon coupled with its rotational motion. The 
integrator has the option of integrating tne plane- 
tary motions or reading their ephemendes from 
DE96* A relatively new integration method is 
utilized in tne UMLTR as developed by E. 

Everhart. This method is an implicit single- 
sequence method using "Gauss -Radau spacxngs for 
substeps within each sequence” . Integration 
orders of 7 to 31 are possible. Other novel 
features of the UMLTR are that Euler parameters 
are used to model the lunar rotational motion and 
the figure-figure interactions (FED of the Earth 
and Moon are modelled from the mutual potential 
function , 

The effect of FFI interactions and coupling 
on tne Moon's motions has had little discussion in 
the literature. A preliminary numerical study of 
the effect of figure-figure interactions on the 
lunar physical Iibrations was given by 3reedlove . 3 
loder 1 * investigated the effect of spin-spin inter- 
actions on the lunar pnysical Iibrations from a 
theoretical standpoint, Lan grange 5 investigated 
the effects of lunar libra tion on the lunar orbit, 
rf. J. Eckert 6 and Spencer Jones 7 attributed dis- 
crepancies between the ooserved and modelled 
motions of the lunar node and perigee to tne 
"lunar figure, inaccuracies m observed and 
theoretical values, and possible effects of un- 
known forces". Based on these discrepancies, 

Eckert estimated values for certain lunar 
-nertias and their ratios. 


In a more general sense, a number of authors 

have discussed tne prooleia of translational/rota- 

tional motion of mutually attracting rigid cocies. 

Dubosnin 3 snowed that if terms factored bv 

* 

r - wnera r is the mutual distance between two 
rigid bodies - are retained in the differential 
aquations then one has the case of semi -indepen- 
dent motions. That is, the translational motion 
can be solved independently and then used as 
input to the rotational solution. If terms 
factored by r are retained then translational/ 
rotational motions are coupled and if terms 
factored by r " 5 are retained the FFI terms appear 
in the torques. Finally, tne FFI terms appear 
m the force expressions if terms factored by 
r " 6 are retained. Beietskn 9 mentions that 
7, T. Kondurar "observed that the perturbations 
introduced in tne lunar motion dv its "shape 
effect" are of the same order of magnitude as some 
of the effects considered m the modern theory of 
the Moon". Several special-case solutions are 
available in the literature cased on assumed mass 
distributions m coth codies . Examples of these 
are given by Lanzano 10 and Kinoshita 11 , 
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More practically/ a vigorous treatment of 
the orbital-rotational coupling may provide for 
improved accuracy in the determination of the 
lunar moments of inertia from observational data. 
This will provide more confidence in the model of 
the lunar interior. 


ORIGINAL PAGE IS 
OP POOR QUALITY 

The purpose of this paper is to present the 
results of some numerical experiments made using 
the UMLTR m an investigation of FFI effects, 
coupling effects and the feasibility of esti- 
mating the parameters of a, 0, y, C^ n and C 
using the LURE data set. ' 0 


The determination of the lunar inertias 
A/MR 2 , B/MR 2 , and C/MR 2 has always been accom- 
plished by combining data obtained from two 
fundamentally different sources. For example, 
the motion of the lunar perigee and node provide 
a measure of K and l wnere 

Ir* (3/2) (C-A)/MR 2 =3C 22 -(3/2)C 20 , 
and 

K=(3-2) (B-A)/MR 2 =6C 22 , 

and the lunar rotational motion provides a means 
of determining 

S=(C-A)/B 

and 

Y=*(B-A)/C. 

Using individually determined L and K values or 
0 and y values leaves tne inertias undetermined. 
But, if e.g., L, 0, K, and y axe algeoraically 
combined, estimates of the inertias can be 
obtained. Michael 12 summarizes various solutions 
for the inertias using the above approach, 

Eckert 5 also utilized this approach. 

Recently, lunar orbiter, Explorer, Luna and 
Apollo orbital data haye been utilized to esti- 
mate 

c 20 = t (A+3) /2-C] /MR 2 

wmch can be algebraically combined with 0 and y 
obtained from tne LURE data m order to solve 
for the lunar inertias, Gapcynski, et al 13 
found C/MR 2 =0 , 392±0 , 003 using this approach, 
Sinclair, et al 14 proposed a simultaneous esti- 
mation of all six lunar inertia tensor components 
using LURE data and Lunar Orbiter IV data. 

Typically, 0 has always been better deter- 
mined tnan C 2Q or Y and Williams 13 reports that 
the uncertainty in C/MR 2 is dominated by the 
error m C 2Q which means that orbiting space- 
craft data must be improved. 

Since the UMLTR models the coupled orbital- 
rotational motions of the Moon, the rotational 
parameters a-(C-A)/B,0, and Y and the orbital 
parameters C 2Q and C 22 =(B-A)/4MR 2 appear in the 
equations of motion. These five parameters are 
not independent however, but are connected by 
the constraints 

C 22“Y U^B) C 2q /2 [Y (1-8) -20] (1) 

and 

a=(0-Y)/(l-3Y> . 

Thus, the use of tne UMLTR in the data reduction 
process offers the possibility of solving for tne 
above parameters simultaneously using a single 
data set sucn as the LURE daca set or tne VLBi 
(Very Long Base Interferometry) data set witn the 
attendant accuracy of such data. 


Section 2 provides a brief summary of the 
UMLTR, Section 3 presents the results of a 
numerical study using the UMLTR which gives some 
indication of the magnitude of orbital -rotational 
coupling effects on the lunar orbit. Section 4 
summarizes previously obtained results concerning 
FFI effects. Section 3 presents the results of a 
preliminary covariance study which indicates that 
the lunar inertias A/MR 2 , B/MR 2 and C/MR 2 can be 
estimated simultaneously using a single accurate 
data set such as that available from the LURE. 

2. The Unified Model of Lunar 
Translation/Rotation- (UMLTR) 


The equations of motion utilized m tne 
UMLTR are based on the general equations. 


m r =V U+F. 
ill i 


(i=l, 11) 


(3) 




o, 1, 2, 3) 


(4) 


wnere is the inertial position vector of the 
ith body (Sun, 9 planets, Moon) , U is the gravi- 
tational potential function, F is a vector of 
forces, (8^} is a set of Euler parameters for 
lunar orientation, and ® } is an angular velocity 


More explicitly, m tne UMLTR, the helio- 
centric orbits of the planets and tne Earth-Moon 
bary center are integrated as subjected to n - 
body effects, relativistic effects modeled using 
the Eddington-Robertson equations as given by 
Anderson 1 ®, and indirect effects due to the 
interaction of the Sun with the Earth figure 

(C 20 and C 30* and MO ° n figtire (c 20' C 22 / S 22 5 “ 

The effect of Ceres, Pallas, Juno and Vesta on 
Mars is also included. 


The geocentric motion of the Moon is inte- 
grated in the UMLTR, as subjected to n - body 
forces, tidal coupling modeled according to 
Oesterwinter and Cohen 17 , Earth figure - Sun 
interaction effects, and relativistic effects 
modeled using the full "N - Body" metric 
suggested by Moyer 18 m the form given by 
Stanaish 1 9 . 

The orientation of the Moon is specified by 
locating its principal axes, {z } , with respect 

to a set of reference axes tZ ) using the Euler 

l 

parameters [0 Q , 0^ 0 2 , 8^J m the form given 
by Morton, et al 20 , Axes izj witn unit vectors 
are defined by K^-i^, *‘2 = ~ 1 \ r K 3 =1 & 
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wnere i 


and i, are unit vectors of a 


r 7 "V 

spherical polar coordinate system locating the 
lunar mass center with respect to a geocentric 
'mean equator and ecliptic of 1950.0 system. 
Matrix [B (6 ]] of equation (4) is an ortho- 
gonal matrix^ of Euler parameters and angular 
velocity vector {ft^} of that equation is the 

angular velocity of the Moon with respect to 

lz > . 

i 


-.9. 7xl0 -4 "/cy 2 and -4.7x10 ^"/cy 2 respectively. 
Also? small bi-monthly terms are present in tne 
node and inclination modulated by the 17 3 -day 
solar nodal passage terms. The eccentricity and 
perigee exhibit a purely periodic character with 
small monthly terms modulated by the 205-day 
solar apsidal passage terms. Finally, a periodic 
growth in se m i -major axis of -5cm. is noted as 
well as maximum growth in longitude residuals of 
.001” and in perigee of «QG07” m 3 years. 


The torques modeled in the UMLTR involve 
those due to particle- figure interactions between 
the Earth, Sun, and lunar second degree harmonics, 
particle-figure interactions between the Earth 
and lunar third and fourth degree harmonics and 
figure-figure interactions between the Earth 

(C 20' C 22 ! and the ! . toon (C 20 ' C 22 ) ’ 

An associated program, EMAM, has been 
developed that provides lunar orientation from 
Eckhardt 1 s 2 1 300 series libration theory and 
lunar position from DE96. 

3 . Translational/Rotational Coupling 

If terms factored by are retained in the 
differential equations, coupling is present and 
the forces depend on lunar orientation. Equiva- 
lently, lunar second degree harmonics appear 
in the forces and lunar third degree harmonics 
appear m the torques. Existing developments, 
e.g., chat of Oesterwmter and Conen 17 contain 
the second degree harmonic terms m the forces. 
Some adopted model for lunar rotation must be 
used in those developments however, to provide 
the lunar orientation. The adopted model could 
range from one based on Cassini's laws to one 
based on analytical or numerical solutions for 
lunar rotation made using an adopted lunar orbit. 
This process is not necessary m the UMLTR how- 
ever, since the lunar orientation is automatically 
provided. 

The effects of coupling are numerically 
investigated using the UMLTR by determining tne 
effect of various lunar gravity models on the 
lunar orbital elements a, e, i, w, ft and true 
longitude, L. Here, eacn run has identical 
translational initial conditions out the lunar 
rotational initial conditions and the lunar 
gravity model differ for each run. Three models 
are used viz,, (i) a lunar second degree harmonic 
model (II) , (n) a lunar second and third har- 
monic model (III) , and Cm) a lunar second, 
third, and fourth degree narmonic model (IV) . 

In each run, rotational initial conditions are 
determined from a least squares fit of the UMLTR 
to Eckhardt's 21 300-301 series for the lunar 
physical librations p, Id, T. These fits are 
made over a 3-year period to supress large free 
librations in the longitudinal libration x that 
appear witn an arbitrary choice of initial con- 
ditions. Translational initial conditions come 
from DE96 and all parameters are consistent witn 
DE96 and Eckhardt’ s 300-301 series. 

Figures 1 and 2 snow tne residuals for a 
3 -year period from the Julian day epoch 
2441200.5 between models III and IV Of prime 
significance are the -1.5 ,T /cy 2 quadratic term 
' in longitude wmcn corresponds to an a=0.27cm./cv. 
The quadratic terms in node and inclination are 


The apparent growth m all of those curves 
are in reality, a portion of a longer period 
effect tnan the 3-year integration time. 

The residuals oetween models II and III 
exhibit the same behavior as shown in Figures 1 
and 2. The magnitude of the effects are much 
larger, however. The largest effects are growths 
m semi -major axis longitude and perigee residuals 
of -7m, -0 17” and -0.1”, respectively, as well as 
a quadratic term -160”/cy 2 in tne longitude. This 
latter value is consistent with an a of -24m/cy 
The node and inclination exhibit quadratic terms 
of -0,1 "/cy 2 and~0.07 TI /cy 2 , respectively. 

It is interesting to note that tne ratios of 
the quadratic terms between the II-III residuals 
and the III-IV residuals is the same as the ratio 
of the constant offset m T (longitude libration) 
predicted for those cases by Eckhardt *s 300 series. 

The model in the UMLTR directly includes the 
effects of second order lunar harmonics in the 
forces governing translational motion. The inclu- 
sion of third and fourth degree terms m the 
torque models can only affect the translational 
motion through the coupling mechanism. That is, 
the presence of these torques affects the lunar 
orientation waich, in turn, affects the forces. 

This analysis snows that even fourth degree 
lunar harmonics can have a measurable effect on 
the lunar orbit contributing at least 5 cm. 
effects in the geocentric lunar radius after 3 
years. For consistency, if fourth degree lunar 
harmonics are included m the lunar torques, 
then third degree lunar harmonics should be 
retained m tne forces. 

4 . Figure-Figure Interactions 

FFI arise in the mutual potential of the 
Earth-Moon system wnen the usual assumption that 
Earth and lunar radii are small compared to their 
mutual distance is discarded 

Breedlove 3 considered the interaction of C^ 0 

and C * 22 * or “^ e Earth witn a triaxial Moon using 

the UMLTR. In that paper, a numerical study of 
the effects of the FFI was presented wnich snowed 
that monthly variations of the order of -0,002" 
were present m tne latitude liorations (p and ICT) 
while long period variations in p, IQ, and T 
should be expected with amplitudes -0 01" and 
periods that are multiples and fractions of the 
lunar nodal period (13,6 years) , These effects 
cause range variations of -10cm, Thus, it was 
suggested that these terms be modeled m any 
analysis of LURE data accurate to -0.01”, 
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SuDsecuentiv, Yoder 4 used an approximate 
theoretical analysis to show that the major FFI 
effects produced the following per turnat ions m 
the latitude iibrations, <$p, and the longitude 
lihrations, 6t* 

<Sp=“0V074 exp C-i<p')~0y0026 exp (-iG) (5} 

+ OV0014 exp C-iCJU^'Jl-OJOOl exp (-1 ( 2A-M> * } ] 
<$T=0Y0067 sin (Q-^M'OVOOl sin 2b>. ( 6 ] 


DE96 at the initial epoch JD2 440468 ,9217100694 (UT) 
and B , B were taken from program EMAN at that 

*0 1 o 

initial epoch. Nominal values for tne other para- 
meters were ' 

8 = 6,31x10"** 

Y - 2.275x10-* 

C = 2,075x10'*, 


In the above, is the node angle describing the 
Earth's spin axis with respect to the ecliptic, 

Z is the lunar mean anomaly, (o=(F-L) is the lunar 
argument of p eric enter, and 1 is the imaginary 
unit. 


5, A Covariance Analysis for the 
Simultaneous Estimation of Lunar 


Y '- and C 20 


The feasibility of simultaneously estimating 
lunar 3, y, and C ^ using only LURE data was 

investigated by forming a normal matrix based on 
the LURE observable, T*, numerical partial deriva- 
tives and actual observation times over the period 
September 1967 to January 1973, This normal ma- 
trix was then inverted to give a covariance/cor- 
relation matrix providing a relative indication of 
the feasibility of estimating the above quantities * 

Nineteen (19) parameters were estimated as 
follows : 

0 - lunar geocentric initial oosition vector, 

p - lunar geocentric initial velocity vector, 

B^ (i=0, 1,2) - lunar Euler parameter initial 

“0 values , 

B (i= 0 , 1 , 2 ) - lunar Euler parameter initial 

0 rate values , 

B = (C-A)/B - lunar inertia ratio, 

Y * (B-A)/C - lunar inertia ratio, 

C 2 o = ^ a+ 3 V 2 -C]/!IR 2 - lunar oblateness, 
a = (C— B)/A - lunar inertia ratio, 

C ^2 * (B-A)/4Jla 2 - lunar tesseral harmonic 

B- - lunar Euler parameter initial value, and 
.0 

8 , - lunar Euler parameter initial rate value, 

0 

The first 15 of the above parameters were treated 
as independent and the remaining four were sub- 
ject to tne following exact constraints 

Si =a=(B-Y)/Cl-3Y), (?) 

S 2 =C 22 =Y(1^6)c 20 /2[Ya-3)-261 , ( 8 ) 

S 3 =S 3 =: 3 1 2 -6 2 2 ', and (9) 

s 4 = B 3 = - ( 1 / 83 ) (3 0 3 Q -2 1 3 r 8 2 8 2 ) . ( 10 ) 

Ini rial values for , 0 _ were taken from 

G o G o 


Following Moyer 

r=C^ i r 

where 
J= 


A - 
and 


(a +a s pw/a +a s ), 

ix x x / \ a s xy 

M 
[Dp J, 


the covariance matrix is 

(ID 

( 12 ) 


x 

D 


In equation (11) , p is the parameter set, T* 
is the LURE observable, x are the independent: 
parameters , s are the constrained parameters , W 
is the data weighting matrix and is the a priori 

covariance matrix for x. 

Following Mulholland 23 , the elements of 
matrix A are computed from 


3t* ff*. p P G 3 p1 sTIM 

^7- [^T 3p x " 3 pJ\ c / 


(13) 


where C is the speed of lignt, p is the 
•geocentric position vector of tne laser trans- 
mission station and is the selenocentric 

position vector of the retro f lector ooth expressed 
in the 1950,0 mean Earth equator and equinox 
system. Unit vector 

«+■ M- -*- \ I-* 

□*■= I Pq^O +{C ,, j/|p G -p +< J , Tacitly assumed 

in the deviation of equation 13 is that the 
Earth transmission station and the lunar 
retroref lector do not move during the -2.5 
seconds roundtrip light travel tune. 

The derivatives Bx/Dp^ were computed as 
follows : 

1. LURE normal point observation times 
'taken from Shelus 24 were converted from 
UTC to ET using tne formulation given oy 
Moyer 22 and JPL smoothed 3IH data for the 
pole position, Ai-UTl ana Al-UTC as 
supplied by Williams 25 , 

^.n _^n — * 

2. The vectors p , p , K , and D were 

G 

calculated at tne LURE observation tunes 
using the UMLTR, and finally, 

3 Numerical part 1 a Is for 3p /Dp and 
o< /3p 1 were generated using the UMLTR. 
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At eacn ooservation time, the matrix J was 
incremented using a data weight corresponding to 
the typical LUBE data accuracy of 13cm. or 
4. 3xl0~ 1 0 sec. as given by Williams 15 . 

The covariance-correlation matrix was formed 
after 528 days, 387 days, and 1212 days. 

There were no significant correlations 
among 0, y, C and the remaining parameters. 

The formal 3candard deviations in 8, y, and C^ Q 
were as shown in Table 1. 


Standard 

Deviation 

- 

Parameter I 

8 

J 

C 20 

After 528 days 
After 887 days 
After 1212 days 

. 29x10 ' l 3 
.37xl0" 14 
. 24x10 ~ lh 

.21x10 * 7 
, 25x10 ' a 
. 86x10 

. 68X10 ” 7 
. 27x10 ^ 
. 16x10 “ 8 


TABLE 1. STANDARD DEVIATIONS FOR 8, y, C^ Q . 

The aoove figures are meant to provide a 
relative indication of tne feasibility of simul- 
taneously estimating 8, Y, and C ^ and hence, 

A/Ma 2 , B/Ma 2 , and C/Ma 2 using the LURE observable. 
Many other parameters should be included for a 
more definitive study. 

The parameters 8 and y are routinely esti- 
mated using LURE data, one of the latest results 
being presented by Williams 15 as 
8 - (631.26 ± 0,3) x 1CT 5 and 
y = (227,37 r 0,6) x 1(T 5 . 

The LURE observable is more sensitive to 8 
and y as shown in Table 1 than to C 2Q but the 

magnitudes indicate that C 2Q should be deter- 
minable from LURE data in a simultaneous estima- 
tion of these parameters using the (JMLTR. 

6, Conclusions 

Coupling , effects have been shown to produce 
residuals of 5cm., .0007”, and 001” after 3 
years in the orbital parameters: semi -major axis, 

perigee and longitude, respectively. These 
effects stem from the inclusion of fourth order 
terms in the lunar potential wnxcn. modify the 
lunar orientation and, therefore, tne forces 
producing the lunar orbit. 

FFI interactions have been shown to produce 
effects as great as 0.01 1 ' m the lunar physical 
librations. 

Modeling, the dynamics of the Moon in a 
coupled translational-rotational sense as 
accomplished m the UMLTR offers the advantages 
of providing rigorously for tnese coupling and 
FFI effects. It also offers the possibility of 
estimating, for the first time, the lunar 
inertias A/MR 2 , B/MR 2 , and C/MR 2 from a single 
extremely accurate data set vis. the LURE data. 
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